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^T) Abstract 

i<— K |rhe members of the scarce category of Of^ supergiants present properties that are intermediate between regular O-stars and Wolf- 

fvq Rayet (WR) stars. Significant similarities between these transitional stars and WN-type objects are now clearly established, at least 

5_^ in the visible and near-infrared domains, pointing to common stellar wind properties. In this study, we report on the first dedicated 

Q^X-ray observations of HD 16691 (04If^) and HD 14947 (05f^), revealing a soft thermal spectrum in agreement with the expected 

■^ X-ray emission from a single O-type star. However, the X-ray luminosity of our targets is slightly lower than expected for single O- 

type stars, suggesting that the particular properties of their stellar wind has also a significant impact on the X-ray emission of these 

objects on the way to the WN category. We argue that the X-ray under-luminosity of HD 16691 and HD 14947 may be interpreted 

as the signature in X-rays of the intermediate stage between O and WR stars, as a consequence of enhanced wind density. 
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*]pl. Introduction 

o 

$_i A few O-type stars in the Galaxy are considered to be tran- 
c/3 sition objects toward the WN-type, on the basis notably of 
^^ similar spectral morphologies. Among these, one finds super- 
giants with significant emission lines in the i r visible and near- 



infrar ed spectrum such as Olf^ stars (IContil Il976t IConti et al. 



K*" l'1995i) . The probable presence of bubbles around some of these 
■^ stars, reminiscent of t he case of several Wol f -Raye t stars, has 
— .^ also been reported bv ICappa & Herbstmeierl (120001) . Some of 
■r-l- these transition objects have been studied at length in the vis- 
.• ible domain. For instance, three northern Of^ supergiants, 
^ namely HD14947, HD15570 and HD16691, have been in- 
^f) tensively monitored to investigate the i r line profile variability 
T-H in the blue domain (iDe Becker et al.L l2009h . The variability 
^ pattern revealed by the observations could be interpreted in 
•'H the framework a large-scale co-rotating structure in the dense 
rN stellar wind of HD 16691 (04If+), responsible for double-line 
i^ shapes of the prominent emission lines, with a variability time- 
" ■ 'scale of 1.98 d. Such a large-scale structure could be attributed 
to the inte raction of the wind plas ma with a significant mag- 
netic field dDe Becker et al.l l2009b . On the other hand, recent 
theoretical modelling allowed to explain such Une shapes on 
the basis of stellar rotation only (IHillier et al.l 120121) . but it is 
still not clear at this stage whether stellar rotation alone could 
explain de variability pattern revealed by observations. Even 
though the variations exhibited by HD 14947 (05f^) are less 
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spectacular, qualitatively its behavior is very similar to that of 
HD 16691. Despite their interesting properties, stars of this cat- 
egory lack dedicated investigations in X-rays. 

The X-ray emission from single O-type stars is mainly be- 
lieved to be generated through strong shocks intrinsic to stel- 
lar _windSj_asa consequence of the line-driving instability 



jFeldmeier et al.L 119971) . On the other hand, the magnetic con- 



finement of the wind leading outflows from both hemispheres 
to collide close to the equatorial plane may also constitute an 
additional so urce of X-rays, producing a hardening of the X- 
ray emission (IBabel & Montmerleill997h . In this context, only 
a few very early-type stars benefitted from a dedicated obser- 
vation with recent X-ray observatories, such as XMM-Newton 
or Chandra. More specifically, peculiarities of the stellar wind 
of very early supergiants revealed in the optical domain suggest 
that the conditions prevailing in the X-ray active region of the 
stellar wind may also deviate from those of more 'regular' O- 
type stars. We therefore selected as targets, among the category 
of Of^ supergiants that have been intensively studied in the op- 
tical, the stars presenting among the most striking features in 
their stellar winds, i.e. HD 16691 and HD 14947. 

The paper is organized as follows. The observations and data 
processing are decribed in Sect.|2] The results of our analysis 
are presented and discussed in Sect.|3] and our concluding re- 
marks are given in Sect.|4] 

2. Observations 

HD 16691 was a target of the XMM-Newton sateUite 



jJansen et al.L 1200 Ih during the 10th Announcement of Op- 



portunity (AOlO), with the proposal ID 0671 10 (PI: M. De 

Aprill7,2013 



Table 1: Observation parameters. The quoted exposure times (expressed in s) 
are the effective values, after rejection of high baclcground level time intervals. 



Target 


JD 

-2450000 


Rev. 


Exposure time 
MOSl MOS2 pn 


HD 16691 
HD 14947 


5794.938 
6313.458 


2142 
2402 


18640 19360 12227 
18013 18080 15012 



Becker), on 21st August, 201 1, and HD 14947 was a AOl 1 tar- 
get of the same observatory on 20th January, 2013 (proposal 
ID 069088, PI: M. De Becker). EPIC insti'uments were oper- 
ated in Full Frame mode, and the medium filter was used to re- 
ject optical light. In both cases, about 15 % of the exposure was 
affected by a rather high background level due to a soft proton 
flare, and that time interval was rejected for our analysis. Data 
were processed using the XMM-Newton Science Analysis Soft- 
ware (SAS) V. 12.0.0 on the basis of the Observation Data Files 
(ODF) provided by the European Space Agency (ESA). Event 
lists were filtered using standard screening criteria (pattern < 12 
for MOS and pattern < 4 for pn). We note that the aim point of 
the exposure was set to the position of the target in order to 
simultaneously obtain RGS spectra, but the low brightness of 
the two stars in X-rays did not allow us to perform any relevant 
analysis with these data. Our discussion will therefore focus 
on EPIC data. The observation parameters are summarized in 
Table[T] 

Spectra of HD 16691 and HD 14947 were extracted using 
a circular region centered on the target (a: 02'' 42'" 52.03'; 
(J:-h56°54'16.47", and a: 02* 36'" 46.99"; 5: +58° 52' 33.12", 
respectively), with a radius of 30". The background spectrum 
was extracted in an annular region surrounding the source re- 
gion, whose surface area was set to be equal to that of the source 
extraction region. We computed the response matrices using the 
dedicated tools available within the SAS software. We finally 
grouped our EPIC spectra to get at least 9 counts per energy 
bin. 



3. Results 



Considering the results detailed by iDe Becker et al.l (120091) 
for their spectroscopic optical campaign, there is no reason a 
priori to consider HD 16691 and HD 14947 as binary stars. We 
note also that we did not find any hint for the presence of a 
companion through high angular resolution techniques in the 
literature. We will therefore adopt the very reasonable hy- 
pothesis that our targets are single stars throughout our anal- 
ysis and our discussion of the results. As a consequence, we 
will not consider any additional source of X-rays originating 
fro m any putative wind-win d interaction re gion - either ther- 
mal dPittard & Parkinil2010l) or non-thermal (IDe Beckej l2007b 
- likely to contribute to the overall X-ray emission. 
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Figure 1: Simultaneous fit of the three EPIC spectra with a 1-T optically thin 
thermal model between 0.3 and 10.0 keV for HD 16691 (fop) and HD 14947 
(bottom). The residuals of the fit are plotted in the lower panels. 



3.1. Spectral analysis 

The X-ray spectrum of HD 16691 is rather soft, as expected 
for presumably single stars. Therefore, our data do not point 
a priori to a hard emission component attributable to mag- 
netic confinement emission. If such a process is at work, it 
should only be very weakly active. The same statement holds 
for the spectrum of HD 14947. We analyzed EPIC spectra of 
our targets using the XSPEC software (v.l2.5.lij. We used 
composite models including absorption and emission compo- 
nents. For the absorption, we reproduced the photo-electric 
absorption by the interstellar medium (IS M) through a wabs 



compo nent, based on solar abundances from I Anders & Ebihara 



( 119821) . with the hydrogen column density (Nh) as unique pa- 
rameter. In the case of H D 1 669 1 , we used the in trinsic (B - V)o 
(=-0.28) color given by Martins & Plez (2006). Taking (B- 
V) = 0.41 1 as given bv .Markova et al.. (2004), we derived E(B - 
V) = 0.69. Inser t ing th is value in the dust-to-gas relation given 
bv lBohhnetaP (Il978h . we derived Nh = 0.40 x 10^^ cm 2. in 



'http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/, see also lAmaua Jl99q) . 



the cas e of HD 14947, the (B - V) value given bv lMarkova et al. 
(120041) is 0.389, and the same considerations lead to a color 
excess E(B-V) = 0.67, converted into Nh =0.38 x lO^^cm"^. 
The ISM hydrogen column density of the wabs component 
was therefore fixed to those values, respectively for both stars, 
throughout our spectral analysis. The absorption of X-ray pho- 
tons by the wind material of the target is take n into account fol- 
lowing the same procedure as developped bv lNaze et al.l (120041) 
and used in the context of other O -type star X-ray investi- 
gations (e.g. IPe Becker et al.1 120041) . This approach consid- 
ers photo-electric absorption by an ionized material with solar 
abundances. In our spectral modelling procedure, that compo- 
nent will be referred to as wind, with a unique column density 
parameters (Nwind) accounting for all material in the line of 
sight (and not only a hydrogen equivalent column as for the 
wabs component used for ISM absorption). 

The X-ray emission from a single O-type star is expected to 
be produced by thermal processes. We therefore favored mod- 
els simulating the emission from an optically thin plasma de- 
fined by a unique temperature, at ionization equilibrium (apec). 
We did not allow the element abundances to deviate from the 
solar values given bv lAnders & Grevessel (il989.) . HD 16691 



7.4 xlO^Lo^ 2.9x10^'' ergs-' (iMartins et all l2005h . we de- 



and HD 14947 are indeed not sufficiently evolved to present 
abundances altered by the central nucleosynthesis at its surface, 
and the nature and quality of the data do not allow to reach such 
levels of refinement. We fitted a 1 -T model to all EPIC spectra 
(wabs*wind*apec), leading to a typical plasma temperature of 
the order of 0.4-0.6 keV (i.e. about 5-7 x 10^ K). The best-fit 
parameters are quoted in Table|2]for all data sets, for both ob- 
jects. Considering the rather low quality of the data, we will 
focus the discussion of the results obtained from the simulta- 
neous fit of all three EPIC spectra. The spectra and the best-fit 
model are plotted in Fig. [T] for both objects. We note that we 
also confronted our spectra to 2-T models, but it did not im- 
prove significantly the results. 



3.2. Discussion 

The X-ray spectrum of HD 16691 is typical of what could be 
expected for a single O-type star, with a soft thermal emission 
spectrum characterized by a temperature of 0.4-0.6 ke V. On the 
basis of the unabsorbed flux (i.e corrected for the ISM absorp- 
tion) quot ed in Tablej2] and assusming a distance of 2.29 kpc 

(Markova et al.L 120041) . we derived an unabsorbed X-ray lumi- 
331 

nosityof8.7xlO^L 

we obtained a Lx,unabs/Lboi ratio of 2.9 x 10 



nosity of 9.8x 10-^' ergs '. Considerin g a bolometric lumi 



-■ 3.4 X 10^^ erg s"' (IMartins et al 



20051), 



Such a value 

appears quite low when c ompared to the expected ratio that 
should be close to 10"^ dSana et al.1 l2006t iNazeL l2009l) for 
regular single O-type stars. The X-ray luminosity is about 
a factor 3.5 lower than expected, which is significant regard- 
ing the rather tight dispersion of recent determinations of the 
Lx-Lboi relation. HD 16691 appears therefore to be signifi- 
cantly under-luminous in X-rays. In the case of HD 14947, 
the X-ray spectrum is also quite typical of that of a sing le O- 
type star Assuming a distance of 2.29 kpc ( Markova etal.l 
20041) . we converted the unabsorbed flux into an X-ray lumi- 
nosity of 1.2 X lO^'^ergs"'. With a bolometric luminosity of 



termine a Lx,unabs/Lboi ratio of 4.2 x 10" . Here again, this value 
is somewhat lower (~ a factor 2.5) than expected on the ba- 
sis of the Lx-Lboi relation. However, one should remember 
that the luminosity ratios derived above are especially sensi- 
tive to the uncertainties on the distanc e to the stars and o n 
the luminosity calibration proposed by I Martins et al.l (120051) . 
To get rid of these sources of uncertainties, we also deter- 
mined flux ratios, i.e. fx,unabs/fboi, for both stars. Follow- 
ing a standard exctinction law, we corrected visual magnitudes 
for the interstellar extinction on the basis of the E(B - V) val- 
ues used to derive the ISM column density (Sect.[3jjl , and 



we app lied the bolometric correction given by Martins & Plez 
(l2006h to derive bolometric magnitudes. After converting 



these magnitudes into bolometric fluxes, we derived flux ra- 
tios of 8.6 X 10 ** and 7.9 x 10 ^ respectively forHD 16691 and 
HD 14947. These values are slightly larger than those derived 
on the basis of the luminosities, but we caution that the flux- 
based approach is very sensitive to the effect of interstellar ex- 
ctinction that is significant in the direction of our targets . In par- 
ticular, the association of these stars with wind-blown bubbles 
JCappa & Herbstmeien 12000.) might indeed cause some devi- 
ations with respect to the standard extinction law used in this 
study. In addition, the latter approach is still sensitive to uncer- 
tainties on the calibration of bolometric corrections and intrin- 



sic colors adopted bv lMartins & Plezl(l2006h . especially consid- 
ering the status of these stars, with winds significantly denser 
than those of more classical O-type supergiants. In summary , 
and considering the rather tight relation derived bv tNazel(l2009h 
(log (fx,unabs/fboi) - -6.45 ±0.51 for O-typc star in the complete 
EPIC bandpass), we report on a slight under-luminosity of these 
two extreme supergiants in X-rays. Considering the rather high 
local column densities (Nwind) derived from our fits, it is tempt- 
ing to attribute it to enhanced absorption of X-rays by the wind 
material. For the sake of clarity, our luminosity and ffux ratios 
are summarized in Table|3] 

These results should be discussed in the framework of 
the evolution stage of our targets. As mentioned in Sect.[T] 
HD 16691 and HD 14947 are considered to be in transition 
between the O-type and the WN-type. This statement is no- 
tably motivated by the strong similarities reported in the vis- 
ible domain (strong an d broa d emis sion lines) between these 
objects and WN stars tContil 119761) . Su ch similarities hav e 
also been reported in the near- infrared by IConti et al.l (Il995h . 
One may find it relevant to discuss the X-ray under-luminosity 
of HD 16691 and HD 14947 in the same context. Whilst the 
Lx,unabs/Lboi IS to somc cxtcnt fairly well established for O- 
type stars, no equivalent relationship exists for Wolf-Rayet star 
As already discussed for instance by Stevens (2005) . the Lx- 
Lboi relation should result from a simultaneous eff'ect of intrin- 
sic emission and local absorption by the stellar wind material, 
whose respective variations from one star to the other seem to 
compensate fairly well to yield a stable linear relation, at least 
in the case of regular O-type stars. However, the situation is 
quite different when WR stellar winds are concerned. It is in- 
deed now clearly established that single Wolf-Rayet stars pro- 
duce stellar winds with sufficient density and metallicity to be- 



Table 2: Results of the fit of EPIC spectra between 0.3 and lO.OkeV adopting a 1-T model. The quoted error bars are given at the 90 ^ 
observed flux (fx,ob.s) and the flux con'ected for the ISM absorption (fx.unabs) are estimated between 0.3 and 10.0 keV. 



confidence level. The 



EPIC-MOSl 



EPIC-MOS2 



EPIC-pn 



3 EPIC instruments 



HD 16691 



logNwiND 

kT (keV) 

N (cm-5) 
Xl (d.o.f.) 



91 Q4.22.07 

61"" 
3.00^12x10-4 

0.844 (29) 



21.8022-0] 
07511° 

o o. 0-51 . 

1.523:34x10-4 
1.416(26) 



21 9222-13 

4.0l25o«xlO- 
0.902 (99) 



21 8921-*^^ 

0.58o4g 

2.714-2^x10-4 
1.011(176) 



Count rate (cts s 
fx,obs (erg cm-2 s 
fx,unabs(ergcm-- 



') 

-') 
s-i) 



0.011+0.001 
6.55x10-14 
1.55x10-13 



0.010 + 0.001 
5.70x10-14 
1.33x10-13 



0.033 + 0.003 
5.83x10-14 
1.69x10-13 



6.04x10-14 
1.56x10-13 



HD 14947 



logNwiND 

kT (keV) 
N (cm-5) 
xl (d.o.f.) 



21 9022-11 
56°-^2 

2.87|<'4f xlO-4 
1.137(22) 



22 0922-25 

o.3i°:g 

16.25^«„«"xl0-4 
1.013(26) 



91 0922.08 

^^■^^21.80 
510.61 

3.9913388 X 10-4 

0.838(112) 



21 9822-07 

^^■^°21.90 
5.7712^03 X 10-4 

0.880(183) 



Count rate (cts s i) 
fx,obs(ergcm-2s-i) 
fx,unabs(ergcm-2s-i) 



0.009 + 0.001 
5.43x10-14 
1.65x10-13 



0.011+0.001 
6.10x10-14 
2.04x10-13 



0.034 + 0.003 
6.50x10-14 
2.04x10-13 



6.15x10-14 
1.92x10-13 



Table 3: Summary of luminosity and flux ratios for both stars, on the basis of 
X-ray fluxes corrected for ISM absorption. 



Target RL=Lx/Lboi logRL Rf = fx/fboi logR, 



HD 16691 
HD 14947 



2.9x10-8 
4.2x10-8 



-7.54 
-7.38 



8.6x10-8 
7.9x10-8 



-7.07 
-7.10 



come quite opaque to soft X-rays dOskinova et al.L 120031) . A 
striking example is the case of the non-detectio n of the WN- 
type star WR40 in X-rays (IGosset et al.l l2005h . emphasizing 
the strong differences between O-type stars and their evolved 
WN counterparts. This critical change in the X-ray activity of 
early-type stars suggests that a transition should occur during 
the intermediate evolution stage. This is what we may be ob- 
serving in the form of the X-ray under-luminosity of these two 
Of^ supergiants. These facts are in agreement with the recent 
discussion by Owocki et al. (2013) on the Lx-Lboi relation for 
O-type stars. These authors indeed pointed out the expecta- 
tion that moderately optically thick winds of extreme O-stars in 
transition to the WN type could cause some deviations with re- 
spect to the canonical Lx-Lboi relation in the direction reported 
in the present study. 

As Of^ stars are rather rare objects, this issue is still a matter 
of small number statistics. At this stage, on can also mention 
the case of the Olf^ star VB4 in the HMl cluster, also likely 



to be slightly under-luminous in X-rays with a log(fx/fboi) 
of about -7.4 (Naze et al.2013, A&A, submitted). The ob- 
servation in X-rays of a sample of transition objects such as 
HD 16691 and HD 14947 is expected to fill a gap in the cen- 
sus of observational data on early-type stars. In particular, the 
change in the stellar wind properties occurring from the O- 
type to the WN-type, and already revealed in other wavebands, 
should be investigated in X-rays as well. The evolution of the 
properties of stellar winds across the evolution stages of early- 
type stars should alter significantly the capability of these winds 
to absorb X-rays after their production in deeper layers, there- 
fore steepening significantly the Lx-Lboi relation leading pro- 
gressively to the complete extinction of intrinsic X-rays in many 
single WR stars. A large number of early-type stars, mainly in 
the transition zone of the evolutionary scheme, should be ob- 
served to investigate this issue. 

4. Conclusions 

We reported on the first dedicated X-ray observtion of the 
Of+ supergiants HD 16691 and HD 14947, known to be in tran- 
sition to the WN category. The spectrum is compatible with 
a thermal emission and rather soft, with a typical plasma tem- 
perature of the order of 5-7 x 10^ K. However, the X-ray lumi- 
nosity of HD 16691 and HD 14947 is quite low as compared to 
the emission level expected for single O-type stars. We tenta- 
tively attribute this slight - although significant - X-ray under- 
luminosity to the somewhat diff'erent properties of the stellar 



wind of this evolved transition object with respect to regular 
O-type stars. The deficit in X-rays measured in this study is 
indeed reminiscent of the case of single WN stars whose stel- 
lar winds are intrinsically significantly more opaque to X-ray 
photons than O-type winds. We argue that the X-ray under- 
luminosity reported here may be considered as an additional 
criteria supporting the idea that HD 16691 and HD 14947 are 
indeed transition objects, in addition to other criteria invoked in 
the optical and in the near-infrared. 
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